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We search for the decays (pK^-y and -S" — > (l>K^'y in a data sample of 228 milhon BB 

pairs collected at the T(4S) resonance with the BABAR detector. We measure the branching fraction 
B{B- ^ 4>K--y) = (3.5 ± 0.6 ± 0.4) X lO"'* and set an upper limit ^ (?!)]?°7) < 2.7 x 10"*^ 

at the 90% confidence level. We also measure the direct CP asymmetry in _B~ ^ (j)K~'y, Acp = 
(—26 ± 14 ± 5)%. The uncertainties are statistical and systematic, respectively. 

PACS numbers: 13.25.Hw 



Measurements of the branching fractions and CP 
asymmetries of 6 ^ S7 decays provide a sensitive probe 
of the standard model (SM), in which these decays are 
forbidden at tree level but allowed through electroweak 
penguin processes. They are sensitive to the possible ef- 
fects of physics beyond the SM manifesting as new virtual 
particles contributing to loops. These additional con- 
tributions to the decay amplitudes could affect branch- 
ing fractions and CP violation [1]. The SM theoretical 
prediction [3| and experimental measurements [3] of the 
b inclusive branching fraction have uncertainties of 

about 10% and are consistent with each other. Although 
exclusive b sj branching fractions are experimentally 
easier to determine than inclusive ones, calculations for 
the exclusive modes are theoretically challenging due to 
large nonperturbative quantum chromodynamic effects. 
The expected direct CP asymmetry between and 
decay rates in the SM is -(0.1 - 1)% while the time- 
dependent CP asymmetry in neutral CP eigenstates such 
as 5° — > 0^s7 should be a few percent [5|. A signif- 
icantly larger CP asymmetry of either type would be a 
sign of new physics. 

There have already been results published for branch- 
ing fraction and/or CP asymmetry measurements in sev- 
eral exclusive modes: B K*j 

B ri{')Kj and various B Kmrj [9|] modes. 
The Belle collaboration has measured B{B^ (j)K~j) = 
(3.4±0.9±0.4) X 10-*^ and B{B" (j)K"-f) < 8.3 x 10"*^ at 
the 90% confidence level using 96 million BB pairs [lo| . 
We present the first BaBAR. measurement of the branch- 
ing fraction for the charged mode B^ (pK^j and 
a search for the neutral mode B'~' (f)K^^ |21|] using 
228 million BB pairs. We also measure for the first 
time the direct CP asymmetry in the charged mode 
Acp = [N{B~)-N{B+)]/[N{B-)+N{B+)], where the 
flavor of the B is determined by the charge of the kaon. 

The data used in this analysis were recorded with the 
BaBAR detector at the PEP-II asymmetric storage rings, 
in which 9.0 GeV electrons collide with 3.1 GeV positrons 
to produce T(4S') mesons. The BaBAR detector is de- 
scribed in detail elsewhere [llj. Most important to this 
analysis are the tracking system composed of the silicon 
vertex tracker (SVT) and drift chamber (DCH) inside a 
1.5 T magnetic field, the ring-imaging detector of inter- 
nally reflected Cherenkov light (DIRC), and the electro- 



magnetic calorimeter (EMC). The tracking system can 
reconstruct a B decay vertex with a resolution of 70 /im 
along the direction of the beam, and has a transverse mo- 
mentum resolution of 0.52% at 500 MeV/c. The DIRC 
provides kaon-pion separation of at least Aa significance 
for momenta up to 3 GeV/ c. The EMC detects photons 
over an energy range from 20 MeV to 9 GeV, with a reso- 
lution of 2.6% at 2.5 GeV. A detailed Monte Carlo (MC) 
simulation of signal and background processes was per- 
formed using the EVTGEN generator [l3| and the GEANT4 
package [l3j . 

We search for B (j)K^ candidates based on charged 
track combinations and the presence of a high-energy 
photon using a kinematic fitter [l3| to reconstruct the 
intermediate mesons and the B. Each decay vertex is re- 
quired to have a probability greater than 0.1%. Can- 
didates for (f) K^K~ are selected from pairs of oppo- 
sitely charged tracks that have been distinguished from 
pions based on a particle identification (PID) likelihood 
selection algorithm that uses AE / dx and Cherenkov light 
measurements. The same PID algorithm is used for the 
single K~ from the B^ in the charged mode. We keep 4> 
candidates with masses within a ±10 MeV/c^ window of 
the nominal mass [l5| . In the neutral mode, pairs of op- 
positely charged tracks are accepted as candidates if 
they have a combined invariant mass within ±10 MeV/c^ 
of the ATg mass and if the ATg flight length is greater than 
three times its uncertainty. We require the combined (j)K 
invariant mass to be less than 3.0 GeV/c^. In the neu- 
tral mode a veto is applied by removing candidates 
with combined (jjK invariant mass within ±10 MeV/c^ 
of the mass. Photon candidates are reconstructed 
from EMC clusters that are not associated with charged 
tracks, are isolated from other clusters, and have the ex- 
pected photon lateral shower shape. We require an en- 
ergy of 1.5 — 2.6 GeV in the e+e" rest frame (CM frame) 
and we veto photon candidates that form a tt" (77) can- 
didate with invariant mass between 115 — 155 MeV/c^ 
(470 — 620 MeV/c^) when combined with another pho- 
ton of energy greater than 50 MeV (250 MeV). 

We identify signal B decays through the distributions 
of two quantities, missing mass and reconstructed mass, 
that peak around the nominal B mass. The missing mass 
IS ?^miss — \/\Pt(as) ~ P-Bp, where Pt(as) is the T{AS) 
four-momentum and ps is the four-momentum of the 



5 



B (j)K^ candidate after a mass constraint on the B is 
applied. The reconstructed mass m^ec is the B candidate 
invariant mass calculated from the reconstructed energy 
and momentum. We require 5.12 < rrimiss < 5.32 GeV/c^ 
and 4.98 < rrirec < 5.48 GeV/c^. To further discriminate 
B decays from continuum e+e^ ~^ <M {<1 — u, d, s, c) 
background we use two topological quantities: the ratio 
of Legendre moments L2/L0 and the cosine of the an- 
gle between the B candidate and the e~ direction in the 
CM frame |cos0^|. We require L2/L0 < 0.55, where 
Li = J2j \Pj\\(^os9*\\ p* is the CM momentum of each 
particle j not used in the B candidate, and 9* is the CM 
angle between the particle's momentum and the thrust 
axis of the B candidate. We also require | cos0^| < 0.9. 

The selection criteria described above are chosen to 
optimize Ns / ^/Ns + Nb in the signal region, where Ns 
and Nb are the MC simulated signal and background 
yields, respectively, and the signal region is defined by 



5.05 < TOrcc < 5.4GeV/c^, 5.27 < mn 
|cos6'^| < 0.8, and L2/L0 < OAi 



< 5.29GeV/c2, 
Signal MC is 



based on inclusive B Xs"f events generated accord- 
ing to the model of Kagan and Neubert using 
TUb = 4.62 GeV/c^ for the effective b quark mass. Only 
the part of the hadronic mass spectrum above the (f)K 
threshold of 1.52 GeV/c^ is used, with Xg forced to de- 
cay to (f)K. This model does not take resonances into 
account. 

After all criteria are applied, the average candidate 
multiplicity in events with at least one candidate are 1.01 
and 1.07 in the neutral and charged modes respectively. 
If multiple B candidates are found in an event, we select 
the best one based on a formed from the value and 
uncertainty of the mass of the 4> candidate and, in the 
neutral mode, the ifg candidate. The remaining back- 
ground comes from continuum combinatorics, nonreso- 
nant B KK+R-j, B (/jA'tt", and B (j)Kri. 

Signal and background yields are extracted from a fit 
to an unbinned extended maximum likelihood function 
defined by 



C{Ns,NB,a) 



N 



n [NsVs{x^) + NBVB{x^-d)] ■ (1) 



Ns and Nb are the number of signal and background 
events respectively, the index i labels each event in the 
data set, and N is the total number of events used in the 
fit. Vs and Vb are products of the one-dimensional signal 
and background probability density functions (PDFs) for 
each of the observables x — {mmiss, "^rcc, L2/ Lq, cos 6'^}. 
The signal shape parameters are fixed in the fit while the 
background parameters a are allowed to vary. In order to 
fit the CP asymmetries of signal and background in the 
charged mode, the number of B^ and B~ events is deter- 
mined separately: N^ = ^(1 =F A''(jp)nj, where j — S or 

_B, Tij and A'jjp are the total yield and CP asymmetry of 



species j, respectively, and the upper (lower) signs corre- 
spond to the positively (negatively) charged B mesons. 
The signal PDFs for Wmiss and mrec are parametrized 

by 



f{x) = exp 



2cr£ p + aL,Rx'^ 



(2) 



where the parameters aL,R and aL,R determine the core 
width and variation of the width on either side of a; = 0, 
X being the difference from the nominal B mass of m,niss 
or mrpc.. The mmiss background PDF is an ARGUS func- 
tion [131 , with the endpoint calculated event-by-event 
from the beam energy. The Wrec background PDF is 
modeled as a 2"*^ degree polynomial. The signal and 
background models for L2/L0 both use a binned PDF 
with eight bins. The cos 9'^ distribution is modeled as 
a 2"^^ degree polynomial in both signal and background; 
true B candidates follow a 1 — cos^ 9g distribution if the 
detector efficiency is fiat in cos 9g . 

To determine the signal PDF parameters we use a high- 
statistics B^ i4r*°(— > K^TT^)j sample. Once deter- 
mined, these parameters are fixed for the fit to B — > (f)Kj 
data. We determine the selection efficiency by perform- 
ing a fit of the yields on signal MC. 

We apply several corrections to the signal yield and effi- 
ciency before determining the branching fractions. Stud- 
ies of simulated events show that the main sources of 
signal-like (peaking) backgrounds are nonresonant B — > 
KK^K~^ events, and B (pKn^ or B <t)Krj, 
where one of the photons from the tt'^ or 77 decay is 
lost and the other is picked up as the signal high-energy 
photon. We estimate the amount oi B ^ KK~^K~j 
contamination by fitting for the yield in mass side- 
band regions defined by 989 < < 1009MeV/c^ and 
1029 < < 1049MeV/c2. By interpolating into the 
signal region, we find and correct for 0.0 ± 1.5 and 5 ± 4 
events for the neutral and charged modes respectively. 
These contributions are subtracted from the event yields 
determined in the fit. From the known branching frac- 
tion [l8| of B ^ (t)K*{-^ Kn^) we correct for a contam- 
ination of 0.27 ± 0.16 neutral and 1.98 ± 0.32 charged 
events. There have been no branching fraction measure- 
ments of -B ^ (pK-K^ or B ^ (j>Kj]. We assume that 
the branching fraction of the first is no more than one- 
third that of -B — > (pK* and that of the latter is no more 
than B — > (f>K* . Based on this we assign an uncertainty 
of 0.5 neutral and 2.9 charged events due to nonreso- 
nant B (j)K {tt'^ / rf) background. To correct for any 
fit bias, we generate 1000 simulated experiments using 
PDFs with separate components for BB and continuum, 
and embedding signal events from the full simulation. 
The background components are generated using shape 
parameters determined from the full MC simulation. We 
correct for a bias of -|-4.1±0.5 events in the charged mode, 
due to correlations among the observables in signal MC 
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events that are not accounted for in the fit. In the neu- 
tral mode we find a bias of —0.06 ± 0.20, so we apply 
no correction but include 0.20 events in the systematic 
uncertainty of the yield. We correct for efficiency dif- 
ferences between data and MC in charged track, single 
photon, and K'^ reconstruction. These multiplicative ef- 
ficiency corrections are 0.956 in the neutral mode and 
0.975 in the charged mode. The corrected efficiencies are 
(15.3 ± 0.8)% in the neutral mode and (21.9 ± 1.6)% in 
the charged mode, where the uncertainties are systematic 
(discussed below). 

The signal yields, efficiencies, branching fractions, and 
charged-mode CP asymmetry are reported in TablelH We 
calculate the central value of the branching fractions by 



B, = 



(3) 



where i labels either the neutral or charged mode, Ng is 
the corrected signal yield, N^g = (228.3 ± 2.5) x 10^ is 
the number of BB pairs recorded, Ei is the corrected effi- 
ciency, and bi is B{(t) K+K~)[^B{K° tt+tt")] in the 
neutral mode and B{(j) K^K^) in the charged mode. 
The world average branching fractions are taken from 
Ref. [H. We measure B{B- (f>K-j) = (3.5 ± 0.6 ± 
0.4) X 10~6 and B{B° (I^K^j) = (1.3±1.0±0.3) x 10"^. 
In the charged mode we measure Acp = (— 26± 14±5)%. 
In Fig. [T] we show fits to the data projected onto mmiss 
and TOrec- In all cases, the displayed distribution is 
created with the signal region selection applied to all 
other fit variables. We determine the consistency of the 
branching fraction measurements with the assumption of 
isospin symmetry using 1000 simulated experiments in 
each mode with the number of signal events determined 
by the average branching fraction, — 2.8 x 10^^. 
From the distribution of the differences in branching frac- 
tion between the modes wc find an 8.9% probability to 
measure a difference greater than or equal to that ob- 
served in data. 

For the neutral mode we compute the 90% confidence 
level upper limit on the branching fraction. We use a 
Bayesian approach with a flat prior probability for the 
branching fraction in the physical region < < 1 and 
zero elsewhere. As the likelihood (Eq. [1]) is a function of 
several parameters, we determine its dependence on Ns 
by fixing Ng to a series of values and recomputing the 
likelihood at each one, allowing Nb and a to be reopti- 
mized to obtain the maximum likelihood at each point. 
We convolve this function with a Gaussian distribution 
of width equal to the systematic uncertainty of the yield. 
Similarly, for the efficiency uncertainty we also use a 
Gaussian distribution of width equal to the efficiency sys- 
tematic uncertainty. We determine the branching frac- 
tion upper bound Bub from the following expression: 

'''' CiB)dB/ [ C{B)dB^ 90%. (4) 
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FIG. 1: Missing mass (a) and reconstructed mass (b) fits in 
the signal region for the charged mode and the neutral mode 
(c,d). The dotted curves show the background contribution 
while the solid curves show the sum of signal and background. 



After applying the previously discussed corrections to the 
yield and efhciency, and including systematic uncertain- 
ties, we obtain B{B° (j>K°-f) < 2.7 x 10"^. 

We assign a systematic uncertainty to the yield due 
to the fixed signal parameters in the fit. We vary these 
parameters within the ranges allowed by the K*"f sample 
to determine the total uncertainty of the yields. We ac- 
count for other systematic uncertainties due to efficiency 
differences between data and MC in charged kaon track- 
ing, kaon PID, and K^, (f>, and photon selection. There 
are small uncertainties assigned to the L2/LQ selection 
and the tt'^ / rj veto, also due to data-MC efficiency dif- 
ferences. 

Figure [5] shows the efficiency-corrected (f>K invariant 
mass distributions, using the background subtraction 
technique described in Ref. In the charged mode, 

we find that no more than 50% of the spectrum in the 
1.6 — 3.0 GeV/c^ range can come from the 1^2(1770) res- 
onance, and we use this information to bound the uncer- 
tainty due to the assumed MC (pK mass spectrum. We 
determine what the efficiency would have been if half of 
the mass spectrum came from resonant ^2(1770) (j>K 
production, while the other half came from the signal 
MC model. We assign the relative efficiency difference 
between this and the nominal model as an uncertainty. 
Adding all of the previously discussed uncertainties in 
quadrature, we find a total multiplicative uncertainty of 
5.2% in the neutral mode and 7.1% in the charged mode. 
The complete systematic uncertainties for each mode are 
summarized in Table HTl 
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TABLE I: Summary of the branching fractions and direct CP asymmetry. In B^^ 
hmit is also given. 



4>K^^ the 90% confidence level upper 



Decay Mode 



B- 



Yield Efficiency (%) B(IO" 



Acp (%) 



85 ± 15 ± 7 21.9 ± 1.6(syst) 3.5 ± 0.6 ± 0.4 -26 ± 14 ± 5 



;±6±2 
< 16 



15.3 ± 0.8(syst) 



1.3 ±1.0 ±0.3 
< 2.7 



u 

V 

1- 50 
d 



(A 

c 

V 

> 

LU 



1 1 


1 1 


(a) 


1 1 


1 1 



TABLE IL Summary of the systematic uncertainties. Except 
where noted, all uncertainties are given as percentages. 



1.6 1.88 2.16 2.44 2.72 3 

(^-K Mass (GeV/c^) 




(|)-K Mass (GeV/c^ 





Uncertainty (%) 


Source 






K K~'y Subtraction 


19.7 


5.2 


Peaking Background 


6.4 


3.4 


Fit Bias 


2.6 


0.6 


Fit PDF Parameters 


+7.0 
-5.9 


+5.9 
-5.2 


Yield Uncertainty 


events 


-e g events 


Kaon Tracking 


2.8 


4.2 


Kl Efficiency 


1.5 





(j) Efficiency 


1.7 


1.7 


Particle ID 


2.8 


4.2 


Single Photon Efficiency 


1.8 


1.8 


Photon Spectrum Model 


0.4 


2.6 


L2/L0 Cut 


1.2 


1.2 


7r"/r? Veto 


1.0 


1.0 


Efficiency Uncertainty 


5.2 


7.1 


BB Counting 


1.1 


1.1 


Total 


+ 23 
-22 


±11 



FIG. 2: The background-subtracted and efficiency-corrected 
(j>K mass distributions (points with uncertainties) for the 
charged mode (a) and the neutral mode (b). The signal MC 
model for the mass spectrum, based on Ref. p^ . is shown as 
a histogram without uncertainties and is normalized to the 
efficiency-corrected signal yield obtained in data. 

For the direct CP asymmetry measurement we bound 
the K~^/K~ efficiency asymmetry of the detector by us- 
ing the measured combinatoric background asymmetry, 
which is consistent with zero within an uncertainty of 
1.8%. To account for uncertainty due to various peaking 
background sources we assume that each source can have 
a CP asymmetry of up to ±58%, which is the root mean 
square width of a flat distribution between —1 and 1. We 
multiply this by the expected fractional contamination in 
the data sample to obtain the systematic uncertainty. For 
B~ (f)K~' {tt'^ /rj) we assign 1.8% uncertainty, while for 
K^K^-f we assign 3.5% uncertainty. For reso- 
nant B (t)K*{-^ Ktt'^) events, the previous BaBAR and 
Belle measurements [l^ show that the CP asymmetry is 



consistent with zero to within 15%. We therefore con- 
sider it to be negligible. As was done with the branching 
fraction measurement, we vary the fixed signal parame- 
ters of the fit to obtain a 2.2% uncertainty for the signal 
CP asymmetry. Adding the uncertainties in quadrature 
we find a total Acp systematic uncertainty of 5%. 

In summary, we have performed the first BaBAR stud- 
ies of i? — > (pK-f decay modes. We measure B{B^ 
(j,K--f) = (3.5 ± 0.6 ± 0.4) X 10"*^, consistent with the 
result from Belle. We have set a limit B{B^ — > ipK'^j) < 
2.7 X 10~^ at the 90% confidence level. Lastly, we have 
made the first measurement of the direct CP asymmetry 
in B- 4>K-T- Acp = (-26 ± 14 ± 5)%. 
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